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Abstract 

Calculations of the axial- vector component to the radiative correction 
for superallowed Fermi O"*" 0^ nuclear beta decay are here modified with 
quenched rather than free-nucleon coupling constants for the axial- vector 
and electromagnetic interactions with nucleons. The result increases the 
deduced value of Vud but does not restore unitarity in the CKM matrix. 



Superallowed Fermi 0+ 0"*" nuclear beta decays [Q provide both the best 
test of the Conserved Vector Current (CVC) hypothesis in weak interactions and, to- 
gether with the muon hfetime, the most accurate value for the up-down quark-mixing 
matrix element of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, Vud- Recent de- 
velopments 0, ^, however, indicate a deterioration in the quality of the CVC test and 
a lowering of the Vud value such that, with standard values 0] of the other elements 
of the CKM matrix, the unitarity test from the sum of the squares of the elements in 
the first row fail to meet unity by more than twice the estimated error. 

Much of this deterioration is a consequence of the recent improvements P, |^, 
in the calculation of the nuclear-structure dependent part of the radiative correction. 
If the CVC hypothesis were correct, then the J-'t values derived by correcting measured 
ft values for the effects of isospin-symmetry breaking and radiative corrections, should 
be the same for all superallowed Fermi transitions in all nuclei. In the 1990 analysis 
jjjj ^Yie ^^™A£ data point had the lowest J^t value; in particular it was lower than 
the ^^O data point. However the revised radiative correction calculation reverses this 
situation, leaving the ^^O data point with the smallest J^t value. Taken with the other 
seven precision data from ^^™A£ to ^^Co the new analysis suggests a Z-dependence 
in the set of J^t values, where Z is the charge number of the daughter nucleus in the 
beta decay. Such a Z-dependence indicates either an electromagnetic correction is 
still not accounted for or that the CVC hypothesis is false. 

Although we are discussing a purely vector interaction between spin 0+ states, 
the axial-vector interaction does play a role in the radiative corrections. An axial- 
vector interaction may flip the nucleon spin and then be followed by an electromag- 
netic interaction that may flip it back again. This axial contribution to the radiative 
correction was considered by Marciano and Sirlin who cast the result into the 
following expression and estimated its value: 

;^[ln(mp/m^) + 2C] = (0.12 ± 0.18)%. (1) 

Here a is the fine structure constant, rrip the proton mass, mA a mass of order 1 GeV 
that provides a short- distance cut-off, and C represents the nonasymptotic long-range 
correction. We write C as 

C = Cbotii + C'ns, (2) 

where CBom refers to the Born graph in which the axial-vector and electromagnetic 
interactions occur on the same nucleon and Cns is a nuclear-structure dependent 
correction in which the interactions occur on different nucleons. 



1 



In the calculations f^, ||, for C, the axial-vector and electromagnetic vertices 
are evaluated with free-nucleon coupling constants. Yet there is ample evidence in 
nuclear physics that coupling constants for spin-flip processes are quenched in the 
nuclear medium [|], |l^. Thus the purpose of this Letter is to repeat the calculations 
of 10] with quenched coupling constants and investigate to what extent this ameliorates 
the deterioration in the CVC test. 

We assume that the axial- vector and electromagnetic vertices can be described 
by on-shell form factors, even though in the diagrams in question they are off-shell. 
Further we will use nonrelativistic reductions and consider the form factors in the 
zero-momentum limit and characterised by coupling constants glk"^ 0). These 
coupling constants can then be equated with well-known coupling constants of nuclear 
physics as deduced from electromagnetic 7-transitions (and magnetic moments) and 
Gamow- Teller /5-decay transitions p. 



We follow the notation of |10, 11 1 and write the magnetic- moment operator as 



/^S = gili^ + gS.^ + 9il[y2, s], (3) 

where [^2,8] represents a spherical harmonic of rank 2, vector coupled to the spin 
operator, S, to form a spherical tensor of rank 1. Here the superscript, /, denotes the 
isospin structure: 1 = being isoscalar and 1 = 1 isovector. The effective coupling 
constants are written as 

^eff = ^ + Sg, (4) 

where g is the free-nucleon value and 6g a nuclear-medium correction. The free- 
nucleon values are: gi°^ = 0.5, gi^^ = 0.5, # = 0.88, gi^^ = 4.706, gi^^ = 0.0 
and g^^ = 0.0. The nuclear-medium correction can also be expressed in terms of a 
quenching factor 

<l = 9cs/9- (5) 



Calculations of the nuclear-medium correction Sg are given in ||T0|, |Tl| for closed-shell- 
plus (or minus)-one nuclei A = 5, 15, 17, 39, 41. They are based on corrections to the 
single-particle wavefunction for these nuclei being evaluated through to second order 
in core polarisation, and on corrections for meson-exchange currents and isobars. Here 



we will use the values from Table 26 of |]IU[ and extrapolate or interpolate for other 



mass values, A. These values are in good agreement with the empirical values deduced 
by Brown and Wildenthal in fits of shell-model calculations to experimental data 
on magnetic moments and Ml 7-transition rates in sd-shell nuclei. 
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Similarly we write the Gamow- Teller /3-decay operator as 



(GT)efr = fi-LA.cffL + 5'A,offCr + 5'pA,cff['i^2, cr], (6) 

and note the traditional use of the Pauli matrix cr rather than S. The free-nucleon 
values are (/la = 0.0, = 1.26 and g^^ = 0.0. Calculations of 5g can be found in 
||T0| , 11]. Here we will take values from Table 27 ^ of |1^ and again extrapolate or 



interpolate as required. Here the results are not in such good accord with the empirical 
values of Brown and Wildenthal obtained in fits of shell-model calculations to 
experimental Gamow- Teller /?-decay rates. For example, the quenching factors 
from are 0.761, 0.737 and 0.727 for A = 26, 34 and 38 respectively showing greater 
quenching than the values we propose to use, as listed in Table 1, and obtained from 



the calculations of [|1^, |rT|. However the use of these stronger quenching factors for 
sd-shell nuclei in the present analysis would not alter significantly the conclusions to 
be drawn here. 

In Table |l|, we list the quenching factors to be used here. For p-shell nuclei, 
they represent linear interpolations between A = 5 and A = 15; for the S(i-shell 
linear interpolations between A = 17 and A = 39; while for the p/-shell they are 
extrapolated from the A = 41 values using a scaling factor of (A/Al)^-^^ applied to 
6g. 

For the case in which the axial-vector and electromagnetic interactions occur 
at the same nucleon, the radiative correction is universal, i. e. the same for all nuclei, 
and has the value H, |^ 

CBorn(free) = 3(7a#/, (7) 

where / is a loop integral. With the replacements g^ —>■ qAdA and gf"* — > qi^^gi^^ the 
universality is now broken, and the contribution from 1-body graphs is written 



CBorn = 3rf^/+(gAgf -l)3rf^ / 

= CBorn(free) + (g^gf - 1) CBorn(free), (8) 

where the second term becomes part of the nuclear-structure dependence of the ra- 
diative correction. For CBorn(free), we use the value 0.881 ± 0.030 0]. 

For the two-nucleon graphs the operator is complicated and comprises 12 terms 
as listed in Table 2 of 0. Different terms originate in different pieces of the electro- 
magnetic couplings. Terms 1 and 2 are proportional to the isoscalar spin coupling, 
"^^^There is one typographical error in Table 27. The entry for SgA for A ^ 40 Od^^^ should read 
-0.255. 
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Table 1: Quenching factors used in the present study 



Electromagnetic Weak 









Qs 




^(1) 






10 


1.042 


0.897 


1.173 


0.927 


0.878 




14 


1.044 


0.873 


1.201 


0.934 


0.858 


A = 


26 


1.023 


0.869 


1.146 


0.877 


0.835 


A = 


34 


1.026 


0.850 


1.155 


0.870 


0.812 


A = 


38 


1.028 


0.840 


1.159 


0.866 


0.801 


A = 


42 


1.010 


0.862 


1.133 


0.866 


0.824 


A = 


46 


1.010 


0.857 


1.137 


0.862 


0.818 


A = 


50 


1.011 


0.853 


1.141 


0.857 


0.812 


A = 


54 


1.011 


0.849 


1.145 


0.854 


0.807 



Table 2: Revised values for the nuclear-structure part of C obtained through the 
introduction of quenching factors 



Unquenched Quenched 







Cns 


Cns 




- l)CBorn(free) 


A = 


10 


-1.67 ±0.20 


-1.35 ±0.16 




-0.19 


A = 


14 


-1.15 ±0.30 


-0.88 ±0.23 




-0.22 


A = 


26 


0.25 ±0.05 


0.20 ±0.04 




-0.24 


A = 


34 


-0.17 ±0.06 


-0.13 ±0.05 




-0.27 


A = 


38 


-0.10 ±0.10 


-0.09 ±0.09 




-0.29 


A = 


42 


0.50 ±0.10 


0.40 ±0.07 




-0.26 


A = 


46 


0.16 ±0.03 


0.14 ±0.03 




-0.26 


A = 


50 


0.16 ±0.03 


0.14 ±0.03 




-0.27 


A = 


54 


0.20 ±0.03 


0.17 ±0.03 




-0.28 
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gi^^; terms 3 and 4 proportional to gi^^; terms 5, 6, 9 and 10 proportional to gl^^; 
and terms 7, 8, 11 and 12 proportional to g^\ The prescription then is to modify 
the 2-body operator by its appropriate electromagnetic quenching factor and by the 
weak quenching factor, ^a- Numerical results are given in Table 0. 

The data base for superallowed /3-transitions produced by the Chalk River group 
III in 1990 is here updated to include four new lifetimes [0 and four new Q-values 
T^, O]. The methodology for handling the data and the theoretical corrections 



remains the same as that used in [T|, |^ except for the introduction of quenching factors 
in the radiative corrections. The corrected /t-values, J^t, for the eight precision data 
cases are fitted by a one-parameter function 



JFt = J^t{0) = constant (9) 

or a two-parameter function 

J^t = Tt{0) [I + aiZ] (10) 

The results are given in Table both with and without the quenching factors. 
In both cases the introduction of quenching factors improves the fit (reduces the 
X^)- In the 1-parameter fit there is a reduction of 2.0s in J-'t{0) and a concomitant 
increase in Vud, but not enough to restore unitarity. In the 2-parameter fit, the 
intercept J^t{0) is essentially unchanged and hence there is little change in Vud, but 
the slope ai is reduced. This indicates that the quenching factors are responsible for 
about 20% of the putative Z-dependence in the current data base. Thus, in spite 
of the considerations given here, the deterioration in the CVC test in the precision 
superallowed Fermi decay data persists. 

The author acknowledges a conversation at the WEIN92 conference with Prof. 
Yu. V. Gaponov, who insisted a calculation such as this should be performed. 
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Table 3: Fitted values of J^t{0) and oi for the eight precision superallowed Fermi 
/9-decay data, and the deduced values of Vud and the unitarity sum 





Unquenched 
Cns 


Quenched 
Cns 


1-parameter fit 






J^t(O) 


3075.0 ± 3.5 s 


3073.0 ± 3.3 s 




1.91 


1.40 


Vud 


0.9733 ± 0.0007 


0.9736 ± 0.0007 


V^d + VI + 


0.9959 ± 0.0016 


0.9965 ± 0.0015 


2-parameter fit 








3067.4 ± 3.5 s 


3067.1 ± 3.3 s 




(1.4 ± 0.5) xlO-^ 


(1.1 ± 0.5) xlO-^ 




1.17 


0.94 


Vud 


0.9745 ± 0.0007 


0.9745 ± 0.0007 


Vud Vus ~l~ Vu}} 


0.9982 ± 0.0016 


0.9983 ± 0.0015 
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